HEAT CONTENT ASYMPTOTICS FOR RIEMANNIAN 
MANIFOLDS WITH ZAREMBA BOUNDARY CONDITIONS 
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Abstract. The existence of a full asymptotic expansion for the heat content 
asymptotics of an operator of Laplace type with classical Zaremba boundary 
conditions on a smooth manifold is established. The first three coefficients 
in this asymptotic expansion are determined in terms of geometric invariants; 
partial information is obtained about the fourth coefficient. 



1. Introduction 

Let (M, g) be a smooth compact m-dimensional manifold with smooth boundary 
dM and let ^ be a smooth vector bundle over M. Let 

D = -{a'nd-d,dj + b''dk + c) 

be a smooth second order operator over M with scalar leading symbol; we adopt the 
Einstein convention and sum over repeated indices. We assume the matrix {a'-'} is 
positive definite and use the inverse matrix g^j to define a Riemannian metric on 
M. 

We can write the operator D invariantly as follows. There is a unique connection 
V on y and a unique endomorphism _B of so that 

D = D(V, E) = -(a^^ V,Vj + E) . 

The connection 1 form and endomorphism E are given in terms of the derivatives 
of the total symbol of D and the ChristofFel symbols F by: 

E = c — a^^ {diLOj + LOiLOj — WfcFij*') . 

The boundary conditions we shall impose are at the heart of the matter. We 
assume given a decomposition dM — C-n U as the union of two closed subman- 
ifolds with common smooth boundary C-jz H Cx) = S. Let denote the covariant 
derivative of (p with respect to the inward unit normal on dM. Let S be an auxiliary 
endomorphism of V\ctz- We take Robin boundary conditions on Cn and Dirichlet 
boundary conditions on arising from the boundary operator: 

(l.b) e(/.:= (0;™ + ^0)|{c^_s}®0|c^. 

We refer to Seeley ^| for a more general formalism; see also related work of 
Avramidi Dowker |3|H|, and Jakobson et al. 11 . 

Let e~*'°^ be the fundamental solution of the heat equation; u = e^*^^0 is then 
characterized by the equations: 

(l.c) {dt + D)u = 0, u{x; 0) = (j){x), Bu{x, i) = for t > 0. 

The equality u(x; 0) = (j){x) is to be taken in the sense where cj) G C°° {M;V) 
is a smooth section to V which gives the initial temperature distribution. Let 
(j)* e C°°{M;V*) be a smooth section to the dual bundle V* which gives the 
specific heat of the manifold. We denote the natural pairing between V and V* by 
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(•, •). Let dx, dx' , and dz be the Riemannian measures on M, on dM, and on S, 
respectively. We define the total heat energy content of the manifold by setting: 

l3{(j),(l)\D,B){t) ( {e~^^''<j),<j)*)dx. 
Jm 

It is worth putting this in a more classical framework in the special case where 
D = A is the Laplacian and 5 = 0. Let W'^''^{M) be the closure of C°°{M) with 
respect to the Sobolev norm 

11011?= / {\V^\' + mdx. 

JM 

Let Wq'q^{AI) be the closure of the set 

{</> e W^'^{M) : supp(0) n Ci5 = 0} . 

Thus, for example, W^'^^{M) = W^'^{M). Let W^o'^(M) = Wq'qj,j{M). For A > 0, 
let 

N{M,Cd,\) =sup(dim£;A) 
where the supremum is taken over all subspaces Ex C Wqq^{M) such that 

||V0||i2(M) < A||0||i2(M), V<^e£;A. 

Then N{M, 0, A) is the spectral counting function for the Neumann Laplacian on 
M, N{M, dM, A) is the spectral counting function for the Dirichlet Laplacian on 
M, and N{M,Cd,^) is the spectral counting function for the Laplacian acting in 
L^(M) with Dirichlet conditions on Cd and Neumann conditions on dM — Co- 
lt is well known, see for example McKean and Singer [l^ , that since M is compact 
and dM is smooth, N{M, 0, A) is finite. By the variational principle, 

N{M, dM, A) < N{M, Cd,X) < N{M, 0, A); 

consequently N{M, Co, A) is finite and counts the number of eigenvalues less than 
A with the Zaremba boundary condition defined by Cb- Let Ai < A2 < ... be the 
eigenvalues counted by N{M, Co, •) and let {4>i] be a corresponding orthonormal 
basis of eigenfunctions in L'^{M). If </> and 0* are smooth, we can express the heat 
content in terms of the Fourier coefficients: 

(l.d) /3(0, cf>*,D, B)[t) = (<^- ^^)LHM) {<P* , • 

i 

By Parseval's identity, the series converges for </> = 0* = 1 and for all ^ > 0: 
D,B){t) = ^e-*^'{(l,0,)L^(M)}' 

i 

i 

It now follows that the series in Equation ljl.d|) converges for all smooth (p and 4>* 
and for all t > 0. 

We return to the general setting. Adopt the notation established above. 

Theorem 1.1. Let cj) G C°°{M;V) and let (j>* G C°^{M;V*). There exists a 
complete asymptotic expansion (3{(j),(f>* , D,B){t) ^ J2ii>a f^ni4'i 4>* j D, B)t'^^^ where 
the (3n are locally computable in terms of integrals over Ad , Ct>, C-jz, and E. 

We use the dual connection V to covariantly differentiate sections of V* . Let 
the dual endomorphism S on V* define the dual boundary condition B for the dual 
operator D on C-n- Near the boundary, we choose a local orthonormal frame {ci} 
for the tangent bundle of M so that e,„ is the inward unit geodesic normal vector 
field; on E, we assume Cm-i is the inward unit normal of E C Cp. Let indices 
a, b range from 1 through m — 1 and index the induced orthonormal frame for the 
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tangent bundle of the boundary; let indices w, v range from 1 through m — 2 and 
index the induced orthonormal frame for E. Let Lab ■= Tabm and L„„ :— r„i,(m_i) 
be the components of the second fundamental forms of dM C M and S C C-p, 
respectively. Let R be the Riemann curvature tensor with the sign convention that 
^1221 ~ +1 for the standard sphere in R^. 

Theorem 1.2. There exist universal constants Ci so that: 

(1) Poi^,^\D,B) = Jj^^{cP,^*)dx. 

(2) /3i(0>*,Z?,6) = -^/^^(0,r)dx'. 

(3) /32(0, ^\D,B) = ~ J^{D^, r)dx + JcJ{(l>;.n + S<P, ^*)}dx' 

+ !c-u ^l^-- ) - )}dx' + Co /s )d2 . 

(4) /33(0>* , Z?, S) - 3^ + Scj), 0*„, + 

— ^Laaifj), <l)*)-m + (( j^iaa^bb " \LabLab + \Ramam)4'i <j)*)}dx' 
+ /s{((ciim-l,m-l + C2-^im + C^Ly^u + CiS)(t>, 0*) 
+ C5(0, + Ce{(l),(t)*)-m}dz. 

(5) W^e W CO - -i, C3 = C5 = cg - -5^. 

(6) M^e have C2 — ^€4 — 

Remark 1.3. Our methods did not yield ci. They also did not permit us to 
complete the computation of {c2,C4}. 

Theorem 11.21 follows if S is empty from results of [3 ^ ; the new feature here is 
the additional integrands over S present in P2 and P3 and the partial information 
we have obtained concerning these terms. 

Here is a brief guide to the paper. In Section |21 we use invariance theory to 
establish Assertions (1-4) of Theorem 11.21 In Section O we use product formulas 
and make a special case computation to show 

(l.e) Co = -i, ^-C3-C5 = 0, and C2 - ^C4 + cg = . 

2 -y/TT 2 

In Section0]we use special cases on the half-plane to complete the proof by showing 

11 
(1-f) C3 — C5 = — and cg 



2V^' ' 20F' " 

The remainder of the paper is devoted to the proof of Theorem II . II using results 
of [niinilTH]; we believe the methods of Seeley [T71[TH| could also be used. We 
switch focus completely at this stage. Instead of working invariantly and globally 
in the context of Riemannian manifolds using methods of invariance theory, we 
work locally in Euclidean space in a system of local coordinates. In Sectional we 
introduce the function spaces which we shall need. In Sectional we discuss various 
'model' problems and in Section Q we formulate a basic theorem on asymptotics. 
We obtain a number of estimates and conclude the proof of Theorem 1 1.1 1 in Section 
IHlby establishing a slightly more general result (see Theorem 18. 2|) . 

2. Universal expressions for the invariants f3^ 

Let denote the additional invariant defined by integration over E. Dimen- 
sional analysis then shows that can be computed by integrating invariants which 
are homogeneous of weight n — 2 in the jets of the symbol; see the discussion in 
12 El on this point where a similar analysis was performed which studied the interior 
invariants and the boundary invariants for Dirichlet or Neumann problems. Thus 
trivially, /3f — Pf — 0, and therefore the first interesting contribution arises at the 
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/?2 level; this must be a constant multiple of ((/>, 0*). Assertions (1-3) of Theorem 
II. 21 now follow. 

To study the form of the additional boundary integral over S appearing in P^, 
we investigate the local geometry near E. Fix a point zq G E. Let z = {zi, Zm-2) 
be local coordinates on E so 

g(9^,9^)(zo) = '5™ and 9^5(9^, 9i)(zo) = for l<u,v,w<m-2. 

By considering the geodesic flow from E C Cp, we introduce coordinates 

(z,yi) exp^{yie™_i(z)} 

so that 2/1 is the signed geodesic distance from E to dM ; Cu corresponds to yi > 
and Ctz corresponds to j/i < 0. The metric then satisfies: 

5(a^,af) = and gidf,df) = l. 

Now use the geodesic flow of dM in M to introduce coordinates 

{z, yi, 2/2) -> expi^^ y^^{y2e^{z, yi)} 

where 2/2 is the geodesic distance to dM. We then have 

5(9^,9|)=g(5f,a|) = and g{dl 81)^1. 

The only non-zero derivatives of the metric at zq are then given by the second 
fundamental forms: 

Luv{zo) := -^d^9{du,d^){zo), and 
Lm-i,m-iizQ) :== -^d^gidi,df){zo) . 

The structure group is the orthogonal group 0{m — 2) and we apply H. Weyl's 
Theorem on the invariants of the orthogonal group. Assertion (2) now follows 
by writing down a basis for the set of invariants of weight 1 and applying the 
symmetry 

where D and B are the dual operator and dual boundary condition on the dual 
bundle V* , respectively. The usual product and addition formulas then show the 
constants are dimension free and universal. This completes the proof of Theorem 
0(4). 



3. Relations among the universal coefficients 

The universal coefhcient cq of Theorem 11.21 (^) can be determined by a special 
case calculation. Let M+ be a compact convex subset of with non-empty interior 
and smooth boundary 9M+. We suppose that 9M+ contains a closed line segment 
A of positive length. Let M_ be the reflection of M+ with respect to the line 
determined by A. We assume that A = dM+ (1 dM_. We set N = M+ U A/_. 

Let A = —di — 3| be the usual flat Laplacian. Let 

Ctz (M± ) := a, Cv {M± ) := dM± ~ A, 
Cn{N) := 0, C-d{N) := dN = C-d{M+) U Cv{M^). 
We take (p ^ (p* = 1 and 5 = to define u± on M± and ujv on N. Let 

, , f°° sinhffTT — 7)5) 

3.a C7:=4/ . . , ^T^^^^ ■ 

Jo smh(7rs) • cosh(7s) 
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The manifold N has two cusps of angle 27r at dA. The results of [3 |S] can be 
used to compute /3„, while Theorem ll.2l (3) can be used to see: 

(3.b) /3Ar(0,</.*,A,6)(t) 
/3M+(0,r,A,S)(i) 



We have by symmetry that ujvl^, A, B) — u±{(f>, 0*, A, S) if x e M± and t > 
and trivially the normal derivative of N vanishes on A. Thus: 

/3iv(0, r , A, 6) (i) = /3Af+ (<^, 0* , A, 6) (t) + /?M_ (0, 0* , A, B) (t) 
= 2(3M^{c^,cj,*,A,B){t). 

Since dN = Cx)(M+) U C-d{M^), we may use equations (|3.b|l and H3.c|l to see that 
4co = 2c(27r). By Equation (j^aji c(27r) = —1. We may therefore conclude 

Co = -5. 

We use warped product formulae to obtain the two relationships between the 
coefficients given in Equation (|l.e() . Let Mi :— [0, 1] x be the cylinder with the 
usual parameters (r,6). Set: 

dsl -.^ dr"^ + dO"^ , := {0,1} X [0,7r], 

:= {0, 1} X 27r], S :- {0, 1} x {0, 
Ai:=-a2-92^ ^i:=0, 
01 1, 01 = 1 • 

Since E is discrete, dz is counting measure. Since all the structures are flat, 
/3„(0i,0*i,Ai,Si) =0 for n>3. 

Let e be a small real parameter. Let M2 '■— Mi x and let 8 be the usual 
periodic parameter on the second circle. Let / — f{r, 9) be a smooth warping 
function and define 

dsl = dsl + e^'f^-^'^dQ^ D2 := Ai - e~^'f(^-'^dl, 
dvoh = e^^drdOdQ . 

We take the warping function to vanish identically near r — 1 and focus attention on 
r = 0. Note that D2 is not self-adjoint. We let B2 = Bi induce the same boundary 
conditions; we must adjust 5*2 appropriately to once again take pure Neumann 
boundary conditions on C-jz x as the connection induced by D2 having a non- 
trivial connection 1 form. We let 02 = 1, but we set 02 = e"^-^^'"'^' to compensate 
for the change in the volume element. Set U2 — ui. We verify that U2 — e^*'^^''^^ 412 
by computing: 

{dt + D2)U2 = {dt + Ai)ui = 0, 
U2(r,0,e;O) = wi(r,0;O) = 1, 
B2U2 = 0. 



dx — \ft 



N 



dx' 



ON 



+ t 



Laadx' + 2c{2n) }+0{t2), 



dN 



I dx-Vi\-^[ dx'\ 



+ t- 



Laadx +2co \ +0{t^). 
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Consequently we may compute: 

/3(1,0;,A2,62)W = / U2{r,0,e;t)<j)*2{r,0,e)e'f'-'''''Urdede 



M2 

= 27r / ui{r,e,Q;t)drd9 ^2Trf3{l,l,Ai,Bi){t), so 

J Ml 

(3.d) (33(1,^2,02, B2) = 2^/33(1, 1, Ai, Si) = . 

First take f{r,9) — f{9) to be independent of the radial parameter near r — Q. 
We use equation Hl.a() to see: 

ujr=uje^0, t^e = -f/e, 

Consequently we have 

(3.e) (3f^{cl)2,r2,D2,B) = -^ [ ifegdOde ^ ^ [ fgdQ . 

Note that C3i„„(02, </'2) = -csfeee'^^ and 05(^2, '?!'2);m-i = -c^feee'"^ . Because 

/33^(<^2,0;,i?2,S) = -e / (c3 + c5)/ede, 
we have the desired relationship 

We now take /(r, 6) = f{r) where 



^ - C3 - C5 = . 



/(O) = 0, 9^/(0) 1, and 9^^/(0) = for fc > 1. 

Since B(j)2 = on Cti, only Cv and S are relevant. We follow the discussion in 
Section 3 of [H to show /3f ^ = 0. We have 

OJr—^^, UJ'r = ^, S2 = ■^S . 

Consequently, we may compute on Ct> that: 

|(0;mm,0*) + ^(j), 4>*^^) = 3-, 

= — Wr^r + LOrTger — — , 
^^Laa{4',4'*):m = ~ ^f" i 

(y2"^^*^"^^^ Q^ab-^ab Q-^amam) (y2 6 ^ 6^^ ' 
^3 ~ >vl2 ^ 12 12 ' 12''^ ~ 

This implies that 

- (3f{l,r2,A2,B2) 

{{{c2Luu + CiS)4>2,(l)*2) + Ce{(t>2,(f>2):rn}dz 



= (-C2 + ^C4 - C6)e VOl(S) . 

This establishes Equation (|l.e|) by showing that 

C2 — 5C4 + ce ~ . 
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4. An example on the half-plane 

We now consider an example on the half-plane, where the classical Zaremba 
boundary value problem that we are considering has a simple spectral resolution. 
The two-dimensional Laplacian is given in polar coordinates by 

I d f d\ 1 d'^ 



A = 



r dr \ dr J dip'^ 

We let (j) — define the Dirichlet component and (p — tt define the Neumann 
boundary component. The spectral resolution is then given by: 

(4.a) ^xA^,r) ^ ^J^smiif[k + l/2])Jk+i/2{Xr), keJNo- 

We can can use this spectral resolution to write down the Fourier decomposition of 
the heat content where we assume a suitable decay of (j) and (f)* at infinity to ensure 
this is well defined: 

poo 

f3{<l>,<l>*,V,B){t) Ae-*^'7fe,A((/.)7fe,A(0*) where 

,._n "'0 



fe=0 
oo 







We first perform the A integration using Uni, Equation 6.633, 
dx xe"*^ Jp{ax)Jp{(3x) = — e ~Ip ( — 
This leads to the following representation of the heat content, 

^ Jo 7o J a Jo 

r^+r'^ ( rr' 

(4.b) sin((^[fc + l/2])sin((p'[fc+l/2])e — Ik+i/2 f ^ 

We will choose suitable angular parts for the localizing functions 0, 0* to ensure that 
the angular integrals can be obtained in closed form. As we will see, an arbitrary 
r-dependence can be dealt with. 

The basis for the forthcoming calculation is the integral representation of the 
Bessel function Ik+1/2 where k is an integer (we refer to ^Hl) Equation 8.431.5 for 
details): 

= - r d9e''°''coii{[k + l/2]e) 
Jo 

Sin([fc + l/2]7r) f°^\ -zcoshT-{k+^)r 



' cos{[k + l/2]9) 
(4.c) (-1)'= / dr e-^™^hr-(/c+i)T 



Remark 4.1. If we had studied pure Dirichlet or pure Neumann boundary con- 
ditions, then the relevant Bessel functions would be indexed by an integer rather 
than by the half integer k + 1/2. The second term in the representation given by 
Equation 14.c|l would be absent in such a case. 

Substituting the identity 
(4.d) exp 1-^:^4-^ + ^ cos4 = cxp i-^-^—p^ + ^(cos^? - 1)1 
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into (j4.b|) and using a saddle point argument, we can verify that the first term in 
(|4.cf) is 'responsible' for producing the volume and boundary contributions. We will 
show that the second term is 'responsible' for the contributions concentrated on S. 
As we are interested in the contributions concentrated on S, which we will denote 
by /3^, we only study the second term in H4.c|l . 

We give another derivation of the identity cq = —1/2 to illustrate the general 
idea behind the calculation. We assume (j) and 0* to have the product form 

(/) = r!i(^)i?i(r), (j)* ^n2i^)R2{r). 

1, i — 1,2, and perform the 



We first study a constant angular part ^i{(p) 
angular integrations. 



dip sin If 



1 



This yields the identity 



drr 



+ 1/2' 

dr'r'Ri{r)R2{r')e-'^ 

°° g-(fe+l/2)r 



k=0 



(fc+ 1/2)2 



Thus only r ~ and r' ~ contribute to the asymptotic small t expansion of the 
heat content. 

We substitute y — r/^/t, y' = r' j\Jt, and expand around r = 0, r' = 0. To 
leading order this produces 



/3^(0,,^*,I?,S)(t) 



t 



i?l(0)i?2(0) / dT^(-l 



-(/c+1/2)t 



Ay / dy'yy'e 



fc=0 



_ y +y yy 



(fc + 1/2)2 



cosh T 



The constant cq is determined by this triple integral and the sum over k. We 
perform the j/'-integral using |1U| . Equation 3.322.2; this involves a complementary 
error function which, together with jJU], Equation 6.286.1, yields 



/3^(0,r,I?,6)(i) 



^i?i(0)i?2(0)£^^ 



fc=0 



(fc + 1/2)2 



(4.e) 



^^g-(fe+l/2)r _ 



sinh^ T 



4t cosh T 
sinh^ T 



The integration of the single terms in 1)4. e|l is not possible as is seen from the r — > 
behaviour. In order to use jJOl) Equation 3.541.1, 



(4.f) 



r dre-''^ sinh"(/3r) = -^T^S f A _ " « + 1 
Jo 2"+i/3 \2P 2' 



with the beta- function B{x,y), we need to introduce a regularizing factor sinh"^ t 
in H4.e|l . integrate the single terms and perform the limit ^ 0+ at the end of the 
calculation. We obtain 



/3^(^,0*,P,6)(i) 

t °° 
--i?i(0)i?2(0)5] 



(fc + 1/2)2 



2 fc + - + U- + - V^' - 
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with the psi-function ipix) — {d/dx)hiT{x). The first term can be sumraed with 
the aid of the Hurwitz zeta function, 



fc=0 



fc + 3/2 TT 



C being the Catalan constant. The summation over the derivatives of the T/i-function 
is performed using Equation 8.363.8, 

We use this to write 

1" 



k=0 



1 



hm V(-l)^-V 

fc=0 

hm yy 

1=0 k=0 

oo m 

,1™+ E E 



k=0 1=0 + 5 + 5] ) 

1 



1 



m=0 Z=0 



1 

+ + 1/4)2+"^ ~ (/ + fc + 3/4)2+'' 
1 1 



= lim > < 

1/^0+ ^ I (m + 1 

m— ^ 



(77^+1/4)2+" (to + 3/4)2+"' 

TO + 1 m + 1 



/4)2+^ (to + 3/4)2+" 

hm (Ch(1 + 1/4) - Ch(1 + 3/4)) + ^Cff (2; 1/4) - ^((2; 3/4) 
[^^0+ 4 4 



This aUows us to conclude 

3S 



^(1]+^Ch (2; 1/4) - (2; 3/4) = tt + 8C + ^. 



-2^i?i(0)i?2(0), 



which gives us another derivation of the result that cq — —1/2. Since we have chosen 
a constant angular dependence, we do not have < (f), (jf >-m terms. Also, given the 
localizing functions are assumed to be C°°(M; V), we need Ri{r) — Ri{—r), which 
implies {d / dr)Ri{r)\r=o — 0; so we do not have < 0, (jf >;„i-i terms either and we 
need not consider other terms in the asymptotic expansion for this example. 

In order to obtain information about constants C5 and cg, we need to study 
nontrivial angular dependences. The constant cg is studied by looking at 

17i((^) = 1, Vl2{(p) ^ fi\nip . 

Assuming R2{r) = rR'2{0) + 0(r2) as r ^ 0, we have that 0*^1 r=o = ^2(0)- 
As in the previous calculation, we start by observing that 



dip sin (f 



sin = (— 1) 



fc+i 



1 



1 



2fc - 1 2fc + 3 



This allows us to write the leading term of as t in the form 
/3^(0,0*,I?,6)(t) 



2t3/= 



i?i(0)i?^(0)£ 



k=0 



1 



1 



(2A:- l)(2fc+ 1) (2A: + l)(2fc + 3) 

dy / dy yy e « 2 co^ir^ 
^0 
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To proceed as before with [TO], Equations 3.322.2 and 6.286.1, we observe 
/ := / dy' [ dyyy"e~^-^-^^ 

OO .2 /"OO 



^0 

2 d 
sinh T dr 



dy y e * / dye 
Jo 



The 2/-integral is evaluated using jJU], Equation 3.322.2, the resuhing y'-integral 
with ^D], Equation 6.286.1. This shows that 

2J^ d 1 / 3 „ , 2 \ 

/=-- -^2i^i l,:^;2;tanh2r , 

smh T dr cosh r V 2 / 

with the hypergeometric function 2Fi{a,b]c;x). For the particular parameters in- 
volved, the hypergeometric function is 



p A 3 \ 2(-i + 7r^) 

which for x — tanh^ r yields the identity 

/ 3 9 \ cosh^ r , 

2-Fi 1, -;2;tanh-^r =2 (-1 + coshr). 

V 2 / sinh" r 



This shows that 

, ,— [1 cosh r 
/ = -4V7r 



sinh'^ r sinh'' t sinh^ t j cosh^ (^) ' 
and consequently we have that 
/3^(<^,0*,I?,S)(t) 

^ -3-i?i(0)i?^(0) g I _ ^^^2^ ^ - ^ ^^^2fc + 3) 

(4.g) X / dr e-i'+i)^ 

Jo cosh (|j 

We evaluate the integral using ^U], Equation 3.541.8. With /i = + 1/2 and with 
the standard notation JOI 



X + 1 



J 7, — n 



X + k 

k=0 



This shows that 



OO 

/ cosVm ''^ " ^ " ^^^^^^ ^ " ^^^^^^ 

We note that 

1 1 1 



(2fc- l)(2/c + 1) (2A: + l)(2fc + 3) (^i - l)/i(/^ + 1) ' 
which allows us to write 

/3^((/>,0*,I?,6)(t) 

r,,3 OO , X 

(4.h) ~ _li?i(0)i?^(0) gj^-—^^--^ + 2[/?(Ai + 1) - /3(/i)] j 
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Our final task is the evaluation of the sum over k. To this end, we note that 

oo ^ 

> -? Tw 5v = 0, and 

oo 

+ 3/2) - fi{k + 1/2)) = -/3(l/2) 



fe=0 

We may then conclude that 



TT 

2' 



r,v, B){t) ^ --^R,{o)R'^{o)ty^ 



This shows, as desired, that 



C6 



In order to determine C5, we choose 

r2i((p) — 1 and ri2(¥') =cos(p. 
We again suppose that R2{r) — rR'2{Q) + O(r^). Then as r ^ 0, 

0*™_i|r=o = ■^2(0)- 
The relevant angular part integration is 



(i(/?sin Lp 



.4 



1 1 

cos if = 



2k -I 2fc + 3' 



The y and y' integration, as well as the resulting r-integration are the same as 
before, and the equation corresponding to 14. h|) for this example is 

fi 

(3''{cj,,<f>\V,B){t) ~ - i?i(0)i?i(0) 
37r2 

(4.i) X £(-1)'= + ^ + Amp + 1) - • 

/c=o LM i Ai+i J 

The remaining sum over k may be performed with the help of the identities 

Ei^- 4(^^ + 4), f:^ = i(4-.), and 
fc=0 ^ 2 ^ fe=0 + 2 ^ 

00 

5^(-l)'=(fc + l/2)(/3(fc + 3/2) - P{k + 1/2)) = . 

fc+O 

We add these relations to conclude that 
This shows that 

1 

C5 



2^/¥■ 

This establishes Equation (jl.tjl and thereby completes the proof of Theorem ll.2l 
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5. Function Spaces 

The analysis in question is local so we shall suppose M is an open domain in M™ 
with compact closure and with smooth boundary dM . For the sake of simplicity, we 
shall assume that the vector bundle in question is trivial; the analysis is similar in 
the bundle valued case. We write the Robin boundary operator in local coordinates 
in the form 

Ru = ^ a'^{x)aidxj + d{x), 

where a = (cti, . . . , (t,„) is the outward unit normal to dM and d is a smooth 
function on dM. 

Let v be the distance to dM and let x' denote a point on dM. We introduce 
coordinates (x', — > a;' + i^cr on a collared neighborhood of the boundary to express 



(5.a) D{x,d^) = -a{x')d'^ + vb{x' ,v)dl + Li{x' ,v,d.^,)d^ + L2{x' ,iy,d^>) . 

In Equation 1)5. a|l . a is a smooth positive function on dM, b is smooth in a neigh- 
borhood of the boundary, and Li and L2 are differential operators in x' of orders 
1, 2, respectively, with smooth coefficients near dM. 

Near S, M is diffcomorphic to E x B+{e) for some £ > 0, where 

B+{e) := {y = (2/1, ya) : 2/2 > 0, + < e^} . 

We may choose 2/2 to be the normal parameter v and use coordinates {z, y) near E. 
With these normalizations, 

= {(z,y) : y2 = 0,yi > 0} and Ck = {(z, y) : ya = 0, ?/i < 0} . 

We can express the operator D as: 

, . D{x,dx)^ -L{z,dy) + A2iz,y,dz) 

+ J2i=i,2iy^Bz2{z,y,dy) + Aii{z,y,d:,)dy^) . 

In this formulation, An and Bi2, A2 are differential operators with smooth coeffi- 
cients of orders 1 and 2, respectively. Furthermore, L can be expressed as 

(5.C) i(^,5y)=E.,,=i,2^'^W%5.., 

where A'^^ are smooth real valued functions on S such that the matrix {A*-'} is 
symmetric positive definite and 

(5.d) i(z,0 < >^cr^lCP for all z e S and C e 



with some positive xq. The boundary operator R can be represented in these 
coordinates as 

(5.e) R{x',dx) = Ro{z,dy) + yiBi{z,yi,dy) + Ai{z,y,dz) , 

where Bi and Ai are differential operators with respect to y and z of orders 1 with 
smooth coefficients, and 

(5.f) i?o(^,a,) = E,=i,2^^'W5..- 

We blowup E and introduce polar coordinates near S by setting: 

{z,p,0) {z,y) = (z,pcos6', psin6'). 

Note that 9 = defines C-d while = tt defines C-n- 

Let C°°(Ms) be the class of smooth functions on C°°{M — E) which extend 
smoothly to the blowup. In what follows we shall suppose that a function p, which 
defined only locally, has been extended smoothly as a C°°{M^) function which is 
positive outside the original neighborhood. 
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Definition 5.1. For k > 0, let £{k) be the set of functions U G C°°([0,oo)) which 
satisfy estimates of the form 

|5^t/(T)| < Cfc,«,'exp(-K'r2) foj. k = 0,1,2,... and k' e (0, k) . 
The optimal constants Ck.n' define semi-norms 

Pk,K'{U) inf Cfc^K' 

giving a Frechet space topology on £{n). We use this topology to define subspaces 
of £{k) of smooth functions on Cp and on C-jz'- 

£{Cv,K) = C°°{Cv,£{>i)) and k) = . 

Introduce the halfspace IR^ = {y = (2/1,2/2) : 2/2 > 0} with polar coordinates 
{p, 6) for p e [0, 00) and 9 G [0, tt]. 

Definition 5.2. Let U he a smooth function on R'^ — {0}. We say that Z-/ G if 

(1) \d'^pdlU{y)\ < Ck,jp''-'' for p < 1 and for all k, j. 

(2) For large values of p, the function U admits an asymptotic expansion 

(5.g) z^(2/)-E^-^|^ + E^f(y2)x(f)| , 

where Vj G C°°([0, tt]), where x is a smooth cutoff function on M+ which 
equals 1 for small r and for large r, and where G £{x) with " + " 
corresponding to 2/1 > and " — " corresponding to yi < 0. 

The asymptotic expansion in Definition 15.21 is to be understood in the following 
sense. For any A'^ = 1,2, . . . and for any multi-indices a and 7 with I7I < \a\, we 
have a constant C, which is independent of 9 and of p, so that: 



(5.h) 



JV-l 



y^d^ I U{p, 9)-J2p-U^-^ + Y^ Uf (2/2)x(f ) 



<Cp-^,p>l 



One verifies that the class given in Definition 15. 21 is independent of the particular x 
chosen. 

Definition 5.3. Let A^(i;) = C°°(I]; A^) be the set of all functions U = U{z,y) 
from C°°(I] X (R^ — {0})) belonging to K'^ for every z G S. The coefficients Vj and 
in the asymptotic expansion | |5.g| ) and in the inequality Ij5.h|) may depend on 
z G S. We assume that these coefficients belong to C°°{Y.,£{k)) and that \b.g\ and 
"hjl can be differentiated with respect to z. 



Definition 5.4. Let be the set of smooth functions on M+ so that: 

(1) |5^'V(p)| < Cfcp^-^- forp< 1, 

(2) V{p) - X;^o O'iP"^ &S. p^oo. 

Remark 5.5. Condition (2) of Definition 15.41 means that we have estimates 

N 

\d%V{p)-Y,a,p-^)\<CN,up~''-''-' 
j=o 

for all N and k and for p > 1. 
We put 

7^^(S) :=C°°(I],i?^). 
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6. Model problems 

We shall first consider boundary value problems on a half-line and then subse- 
quently consider boundary value problems on the half-plane. We begin with the 
Dirichlet problem: 

(6.a) ^ifc/2-if(_^), C/(0) = G, 

where 1/ > 0, t > and fc = 0, 1, . . . The function U — U{v) then satisfies 

U" {v) + 2vU'{v) - 2kU{v) = -AF{v) for v>Q, U{0) = G. 
The homogeneous equation for U (with F ~ Q) has two solutions 

^k{i^) = Sr^^ - vy'er''^ ds and (\)k{v) = e-'^'S^e"' . 

Therefore if F = the only solution to (|6.al) decaying for large v is the function 

C/M = 2&^fc(i.)/r((fc + l)/2) . 

Let G = 0. We impose suitable decay properties on F to ensure the following 
integrals converge and set: 

[ T\ lo {s)F{s)ds + C V'fc {s)F{s)ds) ,k odd , 

[ ^WkTWrnm e'''Ms)F{s)ds , k even . 

Proposition 6.1. Let a E C°°{C-d) he a positive function, let F e £(Ct>,k), let 
G € C°°{Gt>), and let k < imiix' ec-D a,^^{x'). Then there exists U G f (Cd, k) so 

{dt~adl)t^'^U{:^)^t^'^-^F{:^) and {/(O) = G . 

Proof. We make the change of variable r = ai to reduce the problem of Proposition 
IS.ll to that given in Equation (|5.a|l with F — F{x') and G = G(x') for x' £ Gp. One 
can then use the formulae given above to see that U £ £{Cv, k) as claimed. □ 

A similar argument can be given to deal with the Neumann problem: 

Proposition 6.2. Let a £ G°°{C-r) he a positive function, let F £ £{Gti,k), let 
H £ C^{Cti), and let k < minxi^c-R. a.~^{x'). Then there exists U £ E{Cti, k) so 

{dt - adDt'^I^Uix', ^) = t^'l^-^Fix', ^) and U'{0) = H . 

Next we study a model problem in the half-space 

(6.b) (5,-A,)(t'=/2z^(^))-t^-/2-i^(_^) for y e audi >0, 

(6.c) U{yi,0) = g{p) for yi > and pdy,U{yuO) = Hip) for yi < 0, 

where fc is a nonnegative integer, where T £ A^, and where G, Ti. £ R^. Equation 
(I6.bll can be rewritten as 



(6.d) {d^p + ^^dp + jjd^g)U + i2pdp + ~2k)U = -AT onR^. 

We adopt the notation of Definition 15.41 We omit details of the proof of the 
following Theorem as it is analogous to the proof given for Proposition 2 |13| . 

Theorem 6.3. Let n £ (0,1) and p, £ (-1/2, 1/2). Let T £ A^f^^ and g ,n £ R^.. 
Then there exists a unique solution to Equations (|6.c|l and (|6.d|l belonging to AJf . 

Next we consider the case when the operator and the right-hand sides in Equa- 
tions (|6.cl) and 1)6. d|) depend on a parameter. Let i?o be the operator of Equation 
(|5.f|) . let L{z,dy) be the operator of Equation (|5.c|) . and let >€q be the constant of 
Equation IS.dfl . 
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Theorem 6.4. Let k e (0, >^o] and ^ e (-1/2,1/2). If T e A(f-2(x;) anrf if 
G jTi (z i?p(S]), i/ien i/iere exists a unique element U G A^(E) smc/i t/iai 

(at ~ L(z, a,))(t'=/2Z^(z, ^)) = t'=/2-i.F(z, ^) /or zj e and t > 0, 

^^{z,yi,0) = g{z,yi) foryi>0, 

p{Ro{z,dy)U)){z,yi,Q) =n{z,-yi) for y^ < . 

Proof Let B = B,j = A'^/'^ where A = A'K Set Yu = Ej=i 2^fejyj - Then the 
equations given in Theorem 16 . 41 become Equations (|6.b|l and 16.c(l where the right- 
hand side depends on the parameter z e S. The desired result now follows from 
Theorem ESI □ 

7. A Theorem on asymptotics 
We can now establish the result from which Theorem II. II will follow. Let 

K = min x(x') 

x'edM 

where >f(x') is the best constant in the inequality 

<^"'(2^')ICP for eelR". 

Theorem 7.1. Let $ e C°°(Me). Then the weak solution of Equation (|l.c|] has 
the asymptotic representation: 



k=0 



/n f/ie above, ^, anc? are smooth cutoff functions on R+ which vanish when 
the argument is greater than 5 and which are equal to 1 if the argument is less 
than ^5 where 6 > is suitably chosen. We set S 1 — The functions uu 
belong to C°^{M^), Uj^ E £{Cv, k), Uj^ G £{Ctz, k) andUk e A(^(i;) with arbitrary 
jJL G [0, 1/2). Moreover the coefficient Uq equals zero. 

The asymptotic expansion in this theorem is to be understood in the sense that 
the difference 

N 

r^{x,t) = u{x,t)-J2t'^/'{E[^)(u,{x)+pv{i^/p)U^[x\^) 

(7.a) + prji^/p)U^(^x',^))+ap)U,[z,^)} 

satisfies the estimate 

(7.b) \drd^r^ix,t)\ < C| ^(^.+i_,)/2-„„^,H„| ^ ; ^2 

Moreover, the same estimate (in a neighborhood of E) is valid for all the derivatives 
of rjv with respect to z where the constant involved may depend on the number of 
derivatives. 

7.1. The main term. We are looking for u in the form 

u{x, t) = S(^)Mo(a;) + v{x, t), 

where wq = tf> and v satisfies the equations: 

{dt+D{x,d^))v = f iuA/x (o,r), 
(7.c) { v = g onCvX (0,T), 

Rv = h on C-R, X (0,T), 
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with initial condition 

(7.d) w = on M for t ^ 0. 

Here 

(7.e) fix, t) = -E{^)Du„{x) + [dt + D, ^{^^)]uoix), 

where here and elsewhere [•, •] denotes the commutator of two operators. Decompose 
the function u near E in an asymptotic series with respect to p: 

CO 

/(x, t) ~ s (^) /oi (x) + ap) E t-'^'^'^ok ^) > 

where /oi € C°°{M^) and !Fok{z, Y) belongs to A(^(E) and is equal to zero for small 
\Y\. The asymptotic expansion given above means that the remainder 

N 

q^{x,t) = f{x,t) - s(^)/oi(x) -Y^r'+'/'j^ok^z, ^) 
satisfies the estimate 

\d^'>d^d]qNix,t)\ < |Ct(^-l)/2-"o-|a|/2^ 

where qN{x, t) — for p < £\Jt and for some small positive e. 
Analogously, one can represent g and h as 

g(x\t) = s(^^)p.goi(a;') +C(p)2(^).goo(^), 
where 500 G C~(I]) and goi e C°°(Cp), and 

CXD 

where /loi G C°°{Cti) and 7ioA;(-z,F) are smooth functions from R^C^) which are 
equal to for |y| < e. 

Now, the function Uq can be found by solving the problem 

(9, + L(z, dy))llo{z, 5^) = t-i(S(5^)(^)-2«oo(^, + ^oo{z, ^)), 
Z^o(2:,yi,0) = E{^)goQ{z) for j/i > 0, 
p{RoUo)iz,yi,0) =noo{z,^^) for j/i < 0. 

By Theorem l6.4l this boundary value problem has a solution from A^(S). Similarly, 
the function [/(f satisfies the relations 

{dt - ^ 0, C/ff (0) = goi(x') 

and, by Proposition 16. II has a solution [/ G £{Cx>,k)- The remainder 

^/;(x,t) = - S(^)K(x) + v{u/p)pUo^{x', ^)) - CW(^) 

satisfies the Equations 

r (at + Z?(a;,9,))w = / inAfx(0,T), 
(7.f) { w = g onCvx (0,T), 

[ i?w = /i on Ct^ X (0,T), 

and the initial condition 

(7.g) w = on M for t = 0. 
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One can verify that the right-hand sides in (j7.t|) admit the foUowing asymptotic 
expansions as t — s- 0: 

oo 
A;— 1 

(T.h) + pr;(;./p)F,«(x',^))+C(p)^fe(z,^)} 

with fk e C°"(Ms), e SiCv, e f (Ck, and Tk e K'H^), 

with gfc e C^iCv), Gk e 7^^+l(S), and 

Vt 



(7.j) hix\t)^f2 1^'-'^/' (^) ^'^ (^') + C(P) f ^) 



with hk G C°°{C-ji), Ti-k G 7?.^+i(E). The asymptotic expansion of Equation H7.h|) 
is to be understood in the foUowing sense. Let 

N 

/W(x,i)=^t^-/2-l{2(_^)(/,(,)+,,(,/^)^^^^(,',_^) 

then the remainder Rjy = / — /'^^ satisfies the estimate 

r f(jV-l)/2-Qo-|a|/2 ^ < 2 

(7.k) \drdSd2R^ix,t)\<c}^ forJip^ 

for aU ao = 0, 1, . . . and muhi-indiccs a = (ai, •, am) and 7 = (71, . . . , 71,1-2) with 
nonnegative integer components. The derivatives with respect to z are defined 
and should be taken into account only in a neighborhood of E, outside of this 
neighborhood 7 is zero. 

If we denote by g'^'^\x' ,t) and h''^\x' , t) the partial sums in (j7.i|) and ( |7.j| ) from 
1 to and introduce the remainder terms Rgiq = g — g^^^ and RhN ~ h ~ ft,'^^ 
then the asymptotic representations H7.i|l and ( |7.j| ) mean that 

r +(A'+l)/2-ao-|a|/2 c ^ < 2 

(7.1) \drd^,d2R,N{x,t)\ < c| \ ; j;^^ 

and 

r ^N/2~a„^\a\/2 f ^ <• 2 

(7.m) \drd^, d]RM^, i)| < C I ^^^/2_„„_,,/2^M-|«| Z t > p^ 

for all ao = 0, 1, . . . and multi-indices a' = (ai, • • • , ckm-i) and 7 = (71, ... , 7m-2)- 
Roughly speaking the estimates for the remainder terms are the same as the estimate 
for the next terms in the asymptotic expansions Ij7.hl) -( |7j| ). 

The form of the right-hand sides (|7.h|l - frj| are more general than we need but it 
is convenient to consider this more general form in order to unify the construction 
of other terms in the asymptotic expansion for u. 

7.2. Higher order terms. We first describe the construction of the terms ui, C/f , 
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The term ui. We take ui = 2/i. Then 




Decomposing the function ui near E in asymptotic series with respect to p we obtain 
that the right-hand side is asymptoticaUy equal to 

oo 

with /ii e C°°(Me) and Tik € A^(E). Moreover Tik{z,Y) ^ for |y| < e for 
some positive e. So, we have compensated the term containing /i in the right-hand 
side of Ij7.h|) and the discrepancy, which came, can be included in the remaining 
terms in the right-hand side in (|7.hp . We shall denote the new right-hand sides by 
the same letters. 

The term lAi . We find this function from the equation 




supplied with boundary conditions 

Z//i(z,?/i,0) = for yi > 

and 

2 

pY,A=^{z){dy^Ui){z,yi,0)=ni{z,-y^) for < 0. 
The discrepancy in the equation brought by this term is equal to 

(a, + D) (t^'\{p)Uk ^) ) - ^"'/'C(p)^fc ^) 

= t^'\{p){D + L)Uu[z, ^) + t'/^[D, apWu{z, 

Using (j5.b|l and the asymptotic expansion at infinity for functions from the class A 
one can show that the right-hand side of the equation above has asymptotics 

where /u e C^{M^), e £{Cv, n), F« e £{Cn,n) and G Ar'(S). 

The discrepancy in the Dirichlet boundary condition is zero and in the Robin 
boundary condition is 

.|.<--(h(^)..,.', + C(.)|h„(,.^)), 

where hik £ C°°{Ctz), ^ifc G So, one can see that Ui compensates 

the terms JFi, Qi and Hi in the right-hand sides of Ij7.hl) ~( |73t . The discrepancies 
brought by Ui have lower order and can be included in terms in the asymptotic 
expansions (|7.h|l "| [77H ) with k> 2. 
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The terms and . We define the function , x' £ C-p, as a solution of 
the boundary value problem 

(dt - aix')dl)it'/'U^{x',i^)) - t-'/'Ff{x',iy), 

U^{x\Q)^g,{x') 

with x' considered as a parameter. We have 




By Equation H5.a|l . the right-hand side in this equation asymptotically equals 

with Tij G A^^-i(E) and F(] e £{Cv, n). 

The term C//^ is then constructed analogously. Thus, we have constructed the 
approximation of the solution which contains all terms in the asymptotic expansion 
of Theorem 17.11 with k = 1 and which compensates all terms in the asymptotic 
representations Ij7.hp -( |731 ) of the right-hand sides of l\7.i\ with k = 1. Continuing 
this procedure we can compensate terms with k — 2,3, ... . Therefore, if we put 



N 




and — u — u^. The function rjv then satisfies the equations: 

(dt + D{x, d^))rN - fN in M x (0, T) 
TN = 9n on Cp X (0,T), 

RujM = hN on Cn x (0,T), 

with the initial condition = Q on M for t ~Q. The right-hand sides in these re- 
lations admit the asymptotic expansions Ij7.h|) -( |73t , respectively, where summation 
is started from k = N + \. 

8. Estimate of the remainder term in the asymptotics 

In this section, we complete the proof of Theorem l7. ll bv establishing the remain- 
der estimate given in Eauation l7.bl The proof rests on a result fTheorem lS.ll below^ 
obtained by Johansson 12 . First we introduce some additional notation. Let T be 
a positive number and Qt = M x (0, T), = x (0, T) and E^ = Ck x (0, T). 
We introduce also weighted Sobolev spaces. Let £ = 0, 1, . . ., (3 e R. The space 
W'^^'^{Qt) consists of functions on Qt with the finite norm 

ii«ii<.^(Q.)-(/ p'^'-''^ E p'^''W'd:u\^dxdt)"'" 

-'Q'^ \a\<2i 

where we set a := {ao, a) and |a| :— 2aQ + \a\. For s — 1/4, 3/4, . . ., introduce the 
trace spaces W^*'*(r^) with the norm 

\Mw;'-ir?) = iIoM->ml2.f^cy* + Jc^P^''\Hx,-)\^^^^ 
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Here stands for the standard Sobolev space on the interval (0,T). If fc is a 
positive integer, then V^^^^^(Cp) is the space of traces on of functions from 
the space Vp{M) with the norm 

The norm in V'^^^^^(Cp) is defined by 

Analogously, one can define the space W^'*'''(r^). 

The closure of functions from the space W^^'^{Qt) equal to for small t will be de- 
noted by Q (Qt)- Analogously one defines the spaces Vr|o^(r^) and W|j(r^). 

If |a| < 2£ — n/2 — 1 then functions u G W^^'^{Qt) have continuous derivatives 
of order a in Qt and 

This estimate can be obtained from the analogous estimate for functions from non- 
weighted spaces (see 0, Chapter 3) and homogeneity arguments. If w S W^^o^(Qt) 
and |a| -|- 2m < 2£ — n/2 — 1 with a nonnegative m then, clearly, 

(8.a) \drd:u{x,t)\ < . 

The proof of the following result is contained in 12 . 

Theorem 8.1. Let I > 1 be an integer and let (3 satisfy 1/2 < -/3 -f 2^ < 3/2. // 

/ e Wf-'^'-\QT), 9 e wI'^"'-''~"\tJ}) and h € Wl'-^'^''-^'\Tf^) then there 

exists a unique solution u G W'^^q {Qt) to problem (j7.fj) . This solution satisfies the 
estimate 

\\'^\\wf'\QT) ^ c{\\f\\wf--'-'-\QT-) + \\9\\wf-^'^''-^'\T;}) 

(8.b) +\\h\\^2l~3/2,l~3/4.fj.a^. 

Now we are in a position to prove the remainder estimate (|7.bp . According to the 
construction of the terms in the asymptotic expansion given in Theorem 17.11 (see 
the end of Sect. I7.2|l the remainder H7.a() satisfies the boundary value problem (|7.f|l . 
( |7.g| ), where the right-hand sides admit the asymptotic representations (j7.hll -( |731 ) 
with summation starting with k = N + 1. Therefore these right-hand sides are 
estimated by the right-hand sides in l|7.k|l - H7.m|l . This implies that the derivative 
of order k with respect to t and all derivatives with respect to z (in a neighborhood 
of E) belong to 

^/3,o (Qt), i^T) and l^^g (F^), 

for £ < (iV — 2A; — l)/2 and 2(i — (3 < 2 + ^, respectively. We suppose here that ^ 
is an arbitrary number from the interval (0,1/2). Now applying Theorem 18. II we 
obtain that dfr^ together with all derivatives with respect to z (in a neighborhood 
of E) belongs to VF|^/(Qt) for 1/2 < -/3 + 2^ < 1 + m and f < (iV - 2fc - l)/2. 
This implies that in a +ncighborhood of S the integral 

/p''''-''^E|a|<2.P'l"l|9,"5ri9>pdy 

is bounded uniformly with respect to t and z. By the usual imbedding theorem we 
obtain 
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for |q;| <2£-1. Choosing cr = 2£ - 1 - /3 close to 1/2 and then taking fi £ (f, 1/2) 
we can rewrite the above estimate as 

(8.c) \d^d^~-^d]u\ < Ct"'p''-\°'\ 

which is valid for \a\ + 2k + 2m < N — 2 and for arbitrary multi-index 7. 

In order to obtain a remainder estimate outside a neighborhood of E one can use 
.a|) which gives 



(8.d) \d^«d^rN{x,t)\ <Ct"' 

for |a| + 2m < 2£ — n/2 — 1 and for p > e where e is a small positive number. 
Estimates H8.c|l and (|8.d|l imply 

(8.e) |at"°9^rAr(x,t)| < Ci"p'"l-I"l 

for |a| + 2m < N — n/2 — 2 and for arbitrary a G (0, 1/2). 

In order to obtain estimate (|7.b|l for rjv we proceed as follows. We choose an 
integer M > N and represent the remainder term as 

M 

rN=rM+ E t'^'{^{/lj){uk{x)+pv{i^/p)Ue{x',^) 



k=N+l 



(8.f) +p^(^/p)u^(^x\^)) +C{p)U,[z,^)}. 

One can check that all the terms in the summation satisfy estimate (|7.b|l . By 

choosing AI sufficiently large, we obtain estimate (j7.b|) for tm from (|8.e() . The 

proof of Theorem 17. II is complete. □ 

Let Ed and E^ denote the boundaries of Cd and Cr, respectively. Clearly, 
functions from C'^{Ms) may take different values on E^i and E^. The existence of 
the asymptotic series given in Theorem 11.11 is a special case of the following more 
general result: 

Theorem 8.2. Let u he the solution to problem 1) i.cj) and let a G C°°{My,)- Then 
the following asymptotic expansion for u is valid: 

(8.g) / a{x)u{x,t)dx ^^t^ I a{x)uk{x)dx + y2t''^^(tak+t^^^bk+t^/^Ck), 



where 

ttk^ f I Vk{z,9)dzd0, bk^ f w^{x')dx' + f w^{x')dx' 

J Yl J — 7r J Cu J Cji 

and 

Ck^ f hf{z)dz+ I h^{z)dz. 

J Yd J^b. 

Here Uk, Vk, , w^, h^ and are smooth functions whose values at a given 
point depend only on values of ip and its derivatives at this point. 

Proof. We have 

(8.h) / (T{x)E(^^)uk{x)dx = a{x)uk{x)dx- I d^^](j{x)uk{x)dx. 
jm ^l\Jt' Jm Jm ^2^ty 

From a, Uk € C°^{My.) it follows that 

N 

a{x)uk{x) = ^7/fe,(z,0)p^' + 0(p^+i) 
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for each N, which imphes that 



/ ^ ( — ^ ) a{x)uk{x)dx ~ t Cj / / Ukj{z,9)dzd9t^^^ , 
Jm ^2v<^ Jt.J-tt 

where Cj are constants independent of the initial data (f. Therefore, the integrals 
in (|8.h|l give the first sum in the right-hand side of ( |8.g| ) in Theorem 18 . 21 and terms 
of the form tak in | |8.g| l . 

Next, consider the integral 



;.i) 



IM 



(8.j) 



Let us introduce a cutoff function C(2/i) which is equal to 1 for \yi\ < e/2 and for 
1 2/1 1 ^ e, where e is a small positive number. Then we represent (|8.i|) for small t as 

(1 - C(yi))w(^)a(x)!7,^(x', ^)dx 

p p OO p oo 

where p is equal to \/yi + y2 in coordinates y — (yi,?/2)- One can check directly 
that the first two integrals in | |8.j| l have asymptotics 

oo „ 
Y,t('+^^^' quA^')dx' 

and the third integral has an expansion of the form 

oo ^ 

Making change of variables yi = and j/2 = 2^/tY2 we can rewrite the last 

integral in (jSjl as 



4i^/' / / / ^{p)pv{y2/pMz,ViY)U^{z,2VtYi,Y2)dYidY2dz. 
is Jo Jo 

Since £ £{Cd,k), the last integral has the asymptotics 

oo ^ 

j=0 -^^o 

(2) 

where hj^J are integrals with respect to Y2 of linear combinations of functions 
dy-^Uj'^ {z,Yi,Y2)\yi=o multiplied by explicit weights. The term 

is considered analogously. 

It remains to obtain an asymptotic expansion of the term 



(8.k) 



p pOC pOQ 

is Jo Jo ^^^^'^^p^^*' ^) '^^'^y^'^y^ 



Using the asymptotic expansion for the function Uk for large second argument: 
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we obtain integrals similar to the ones just considered. Moreover, one can show that 
the coefBcient vqq is equal to zero because of vanishing of the analogous coefficient in 
the asymptotics of the right-hand side in the equation for the function Uq . Reasoning 
as above we arrive at the required asymptotic representation for these integrals. 
This completes the proof of Theorem 18.21 and thereby of Theorem ll.il □ 
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